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We report on the results of the Chandra ACIS-I observation on the central region 
of the Monoceros R2 cloud (Mon R2), a high- mass star- forming region (SFR) at a 
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distance of 830 pc. With a deep exposure of ~100 ks, we detected 368 X-ray sources, 
■ ~80% of which were identified with the near-infrared (NIR) counterparts. We sys- 

Oh 1 tematically analyzed the spectra and time variability of most of the X-ray emitting 

sources and provided a comprehensive X-ray source catalog for the first time. Using 
"55 ! the J-, H-, and A-band magnitudes of the NIR counterparts, we estimated the evolu- 

tionary phase (classical T Tauri stars and weak-lined T Tauri stars) and the mass of 
the X-ray emitting sources, and analyzed the X-ray properties as a function of the age 
and mass. We found a marginal hint that classical T Tauri stars have a slightly higher 
temperature (2.4 keV) than that of weak-lined T Tauri stars (2.0 keV). A significant 
fraction of the high- and intermediate-mass sources have a time variability and high 
plasma temperatures (2.7 keV) similar to those of low-mass sources (2.0 keV). We 
performed the same analysis for other SFRs, the Orion Nebula Cluster and Orion 
Molecular Cloud-2/3, and obtained similar results to Mon R2. This supports the 
earlier results of this observation obtained by Kohno et al. (2002, ApJ, 567, 423) and 
Preibisch et al. (2002, A&A, 392, 945) that high- and intermediate-mass young stellar 
objects emit X-rays via magnetic activity. We also found a significant difference in 
the spatial distribution between X-ray and NIR sources. 

Key words: ISM: stars — ISM: individual (Monoceros R2) — stars: formation 
— X-rays: stars 

1 



1. Introduction 

Low-mass young stellar objects (YSOs) are classified along their evolutionary phases 
into classes 0, I, II, and III based on their infrared (IR) to sub-millimeter spectral energy dis- 
tributions (Lada 1991; Andre et al. 1993; Andre, Montmerle 1994). Class and I sources 
correspond to the youngest and evolved protostars, while class II and III sources roughly corre- 
spond to classical T Tauri stars (CTTSs) and weak-lined T Tauri stars (WTTSs), respectively. 
Low-mass TTSs (=CTTSs and WTTSs) were first discovered to emit significant X-rays with 
the Einstein observatory (Feigelson, DeCampli 1981; Feigelson, Kriss 1981; Montmerle et al. 
1983). The following observations with ASCA and ROSAT found X-rays from class I sources 
(Koyama et al. 1994; Casanova et al. 1995; Koyama et al. 1996; Kamata et al. 1997; Grosso 
et al. 1997, 2000; Tsuboi 1999; Tsuboi et al. 2000; Grosso 2001). In addition, Tsuboi et al. 
(2001) reported on the hard X-ray emission from class candidates using the Chandra X-ray 
Observatory (CXO). Thus, it is now widely known that low-mass YSOs commonly emit X-rays 
during most of their evolutionary phases. Although it is not clear for class sources, X-ray 
emission from most YSOs is explained by a thermal plasma with temperatures of 1 to a few keV 
and rapid time variances, which are similar to those of the solar X-rays. Accordingly, the X-ray 
origin is generally thought to be solar-type magnetic activity: amplification of the magnetic 
field followed by magnetic reconnection associated with the stellar convection layer and differ- 
ential rotation (Feigelson, DeCampli 1981; Montmerle et al. 1983). However, whether or not 
the X-ray emission differs systematically along the evolutionary phases and stellar masses is 
not yet known. 

On the other hand, the X-ray emission of high-mass main-sequence (MS) stars is de- 
scribed by a thin thermal plasma model with relatively soft temperatures of less than 1 keV, 
and shows moderate variability. Lucy and White (1980) proposed that these features are due to 
the plasma heating attributable to the collision of the accelerated stellar wind with circumstel- 
lar matter. However, our understanding of X-ray emission from high- and intermediate-mass 
YSOs is largely behind the low-mass sources and high-mass MS stars, because the former evolve 
very quickly and the population is small; hence, they constitute fewer or more distant samples 
than those of low-mass YSOs. Since YSOs are generally located in dense cloud cores, not 
only optical/near-infrared (NIR) light but also soft X-rays suffer large extinction. Accordingly, 
hard X-ray imaging instruments are essential to resolve embedded YSOs in the cores. The 
CXO (Weisskopf et al. 2002), with a high spatial resolution, large field of view (FOV), wide 
energy range and moderate energy resolution, provides the best instruments for observing the 
star- forming regions (SFRs). 

The Monoceros R2 cloud (hereafter Mon R2) is one of the closest high-mass SFRs, which 
is located at a distance of 830 pc (Racine 1968). The cloud age is ~1 Myr (Carpenter et al. 
1997), providing numerous YSOs over a wide range of ages from protostars to WTTSs. The 
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core mass is ~1000 M Q (Tafalla et al. 1997), which is two orders of magnitude larger than 
those of the cores of low- to intermediate-mass SFRs. Thus, Mon R2 is a good candidate for a 
systematic study of the X-ray emission along the evolutionary phases and a study of high-mass 
star formation. 

The first X-ray observation on Mon R2 was made with the ROSAT satellite (Gregorio- 
Hetem et al. 1998). However, the spatial resolution was limited to resolve individual YSOs 
in the Mon R2 core. Also, a short exposure (~9500 s) observation and poor sensitivity in the 
hard energy band limited the possibility to detect most of the embedded YSOs. Then, the 
ASCA observation on Mon R2 (Hamaguchi et al. 2000) also could not resolve individual YSOs 
due to the limited spatial resolution of ~ 1'. Recently, Kohno et al. (2002) performed a CXO 
observation and resolved X-rays from many YSOs for the first time. The earlier results have 
already been reported by Kohno et al. (2002) and Preibisch et al. (2002), but concentrated 
only on the central 3'.2 x 3'.2 region. This paper reports on a follow-up analysis of the CXO 
data in the entire cloud core and the NIR cluster of Mon R2. The purposes of this paper are 
(1) to make a comprehensive X-ray catalog, (2) to discuss the X-ray properties as a function 
of the evolutionary phases and the masses of YSOs, and (3) to discuss the characteristics of 
X-ray sources by comparing the spatial distribution of the X-ray and NIR sources. 

2. Observation 

Mon R2 was observed with the Advanced CCD Imaging Spectrometer (ACIS) on board 
the CXO on 2000 December 2-4. The ACIS-I array, comprising four abutted front-illuminated 
X-ray CCDs, covered the entire cloud core and the NIR cluster of Mon R2. The telescope optical 
axis position on the ACIS-I array was R.A. = 06 h 07 m 50f86, Decl. = -06°22'50"0 (J2000). A 
high spatial resolution of ~ 0."5, a large field of view (FOV) of 17'4 x 17'A, a wide energy range 
of 0.2-10.0 keV with a moderate energy resolution (AE ~200 eV at E — 6 keV), and a time 
resolution of 3.2 s provide us with high-quality and uniform samples of YSOs in the cloud. 
We used the level 2 data provided by the pipeline processing at the Chandra X-ray Center. 
Besides, we removed those events due to charged particles and hot and flickering pixels. The 
effective exposure time was then ~100 ks. 

3. Analysis 

3.1. Source Detection 

Figure 1 shows an ACIS-I image of the entire cloud core and the NIR cluster of Mon 
R2 overlaid on a 13 CO (J=l— 0) intensity map (Loren 1977). For source detection, we used 
the wavdetect program in the CIAO package version 2.2 1 in the 0.5-10.0 keV band with a 
threshold significance of 10~ 6 and wavelet scales from 1 to 16 pixels in multiples of \/2. In 

1 See (http://asc.harvard.edu/ciao/). 
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order to study the spectra and timing with reasonable statistics, we discarded the sources with 
X-ray counts of less than three-times the background. We then detected 368 sources in the 
ACIS-I FOV. The X-ray photons of each source were accumulated from a circular region with 
a radius of 0."7-15."4, depending on the point spread function (PSF), which is a function of the 
angular distance from the telescope optical axis. The background photons were extracted from 
a source-free circle of ~ 42" in radius. 

3.2. Correlation with the NIR Sources 

To search for the NIR counterparts of the CXO sources, we used the Point Source Catalog 
in the 2MASS Second Incremental Data Release 2 , and the results of SQIID and IRCAM3 
observations (Carpenter et al. 1997). SQIID was mounted at the Cassegrain focus of the 1.3 m 
telescope in the Kit Peak National Observatory, which simultaneously takes J-, H-, and K- 
band images using 256x256 PtSi arrays (Ellis et al. 1993). The pixel scale is l!'36 providing 
a FOV of 5!8 x 5f8. The mosaic consists of 3x3 frames, and hence the SQIID observation 
obtained a ~15'xl5' mosaic image with the photometric mode centered on Mon R2. IRCAM3 
at UKIRT (United Kingdom Infra- Red Telescope) utilizes a 256x256 InSb array. The pixel 
scale is 0''286, providing a FOV of 1'2 x Y.2. The observation consisted of 3x3 grid mosaics; 
hence, the IRCAM3 observation presented J, H, K, and narrow-band L imaging photometry 
of the central ~ 3'2 x 3'.2 region of Mon R2 more deeply than the SQIID observation. The 
5 cr completeness limits in the A'-band images were 14.5 mag and 16.6 mag for the SQIID and 
IRCAM3 observations, respectively (Carpenter et al. 1997). On the other hand, 2MASS has the 
K s -band limit of 14.3 mag. Because the position accuracy of 2MASS (the r.m.s. uncertainty 
of ~ 0"1) is better than those of SQIID (0"3) and IRCAM3 (~ 0"13), we used the 2 MASS 
frame as the standard frame. The distribution of the positional differences between SQIID, 
IRCAM3, and 2MASS could be approximated by a two-dimensional Gaussian in the AR.A.- 
ADecl. plane with the center at AR.A. ~ 0"65 and ADecl. ~ 0"06. We hence fine-tuned the 
SQIID and IRCAM3 frames so that the center of the Gaussian would come to the zero position 
of the 2MASS frame. 

We picked up the closest NIR source from each X-ray source and did vice versa. Using 
the thus-selected closest pairs of NIR and X-ray sources, we fine-tuned the ACIS-I frame by 
shifting — 0"87 in R.A. and — 0"75 in Deck to match the 2MASS frame. We then searched for 
the closest pairs of NIR and X-ray sources again. Among these pairs, we identified counterpart 
pairs under the criterion that the angular distance between the X-ray and the NIR source is 
smaller than the root mean square of the position uncertainty derived by the wavdetect program 
and the radius of the PSF. 

Finally, 245 X-ray sources were identified with the NIR counterparts in the SQIID and 
IRCAM3 databases, and outside of the FOV of SQIID, additional 45 X-ray sources were identi- 

2 See (http://www.ipac.caltech.edu/2mass/). 
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fied with the 2MASS sources (see figure 2). Assuming a random distribution of the NIR sources, 
the chance coincidence that any of the NIR sources fall in an X-ray source region is estimated 
to be ~0.6%; hence, the false identification is expected to be less than two sources. The fine- 
tuned source positions, background-subtracted photon counts in the 0.5-10.0 keV band, and 
the name of the NIR counterparts are given in table 1. 

3.3. Spectral and Timing Analyses 

For bright 126 sources having more than 50 photons after background subtraction (here- 
after "bright sources"), we performed spectral and timing analyses using XSPEC version 10.0 
in XANADU. 3 The X-ray spectra were fitted by a thin-thermal plasma model with three free 
parameters: the absorption column density (Nn), plasma temperature (kT), and normaliza- 
tion. Since the statistics are limited for most sources, the chemical abundance was fixed to be 
0.3 solar according to previous reports (Kohno et al. 2002; Feigelson et al. 2002; Getman et al. 
2002). The best-fit parameters and relevant physical values are listed in table 2 with their 90% 
errors in parentheses. The distributions of Nn, kT, and Lx of bright sources are shown with 
their mean value and that weighed with the errors (Bevington, Robinson 1992) in figures 3, 4, 
and 5, respectively. 

We then made light curves of the bright sources in the 0.5-10.0 keV band. To study the 
time variability, we applied the Kolmogorov-Smirnov test (Press et al. 1988); then, 49 out of 
126 sources were found to be time variables with more than 99% confidence. The significance 
levels of the time variability are listed in table 2. Most of the variable sources exhibit a flare-like 
event with a fast rise and slow decay. 

4. Discussion 

4-1. Cloud Members 

We examined whether or not the detected 368 X-ray sources are cloud members by 
checking their values and i^-band magnitudes. The absorption column density of the 
interstellar medium to Mon R2 is estimated to be ~2.6xl0 21 cm -2 , assuming the density of the 
interstellar medium to be 1 cm -3 . Since the values of most X-ray sources exceed this value 
(figure 3), they should be located in or behind the cloud, although the limited photon counts 
can not give firm constraints. Most of the background galaxies have K > 17 mag (Tsujimoto 
et al. 2002), while the i^-band limit of NIR sources is 16.6 mag. We thus suspect that most of 
the NIR identified X-ray sources are cloud members. 

For all of the ACIS sources, we searched for possible binary systems by checking whether 
or not two or more NIR sources were located within a radius of three-times the PSF. We used 
2MASS and IRCAM3 catalogs, which have comparable pixel scales of 1"0 and 0''3 compared 

3 See (http://heasarc.gsfc.nasa.gov/docs/software/lheasoft/). 
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with that of ACIS of 0'.'5, and found that only one source, #151, has two NIR counterparts 
separated by I'.'l (or ~900 AU). This separation, however, is one order of magnitude larger 
than the typical separation of binaries (Giannuzzi 1989); hence, these two NIRs may not be a 
binary system. 

4-2. Comparison with Previous X-ray Observations and Studies 

In order to make a comprehensive discussion about X-ray emitting sources in Mon R2, 
we compared our results with those of previous studies. 

In a ROSAT/PSPC observation, Gregorio-Hetem et al. (1998) found 23 X-ray point 
sources in the FOV of ACIS-I. We searched for ACIS counterparts within a radius of 5'.'0 from 
each PSPC source (position accuracy ~10") and found seven counterparts. We also found five 
additional possible counterparts relaxing the criterion radius to the position error of each PSPC 
source (table 1). For the remaining 11 sources, we did not find significant X-rays. This may 
have been due to a long-time variability, possibly several years, as can be seen in the Sun. 

There are some sources that are missing in table 1, but are listed in Kohno et al. (2002). 
This discrepancy arises from the difference in the source detection algorithm. In Kohno et al. 
(2002), after executing the wavdetect program, they detected 12 additional sources with a visual 
inspection above the 5 a confidence level. However, in this paper, we did not perform such 
an inspection, but excluded 34 sources under the criterion that they should have the source- 
significance level of more than three. Moreover, since we ran the wavdetect only for the total 
energy band image, some sources detected only in the soft band (0.5-2.0 keV) or hard band 
(2.0-10.0 keV) in Kohno et al. (2002) were missed in this study. We also compared our results 
of four intermediate- and high-mass YSOs (a s , IRS1SW, IRS2, and IRS3NE) with those in 
Kohno et al. (2002). We do not report on X-rays from the position of a s , although Kohno et al. 
(2002) found a counterpart within a 1" radius. We certainly found an X-ray source at a distance 
of 0"54 from a s , which is larger than our criterion radius of identification at the position of a s 
of Of'5 based on the error radius derived by the wavdetect program. We, hence, do not regard 
this X-ray source as a counterpart of a s . For IRS1SW, IRS2, and IRS3NE, we obtained almost 
the same results as that of Kohno et al. (2002). Based on the positions in the color-color 
and color-magnitude diagrams (subsections 4.3 and 4.4), IRS1SW is classified as a high- or 
intermediate-mass CTTS, and IRS2 has a larger iT-band excess than CTTSs. Because IRS3NE 
is saturated at the K band, the reliable .fT-band magnitude can not be obtained. However, 
based on the J- and if-band magnitudes, this is classified as a high- or intermediate-mass 
source. 

Preibisch et al. (2002) resolved a multiple IRS3 source into three individual IR sources 
(IRS3A=IRS3S, IRS3C=IRS3NE, and IRS3E=CMDSH278) by their deep observation using 
HST/NICMOS, and found that all of the sources emit X-rays using the same dataset as ours. 
We found five photons at the NIR position of IRS3S. However, these photons may not be 
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from IRS3 because of severe contamination from the bright source #207 (=IRS3NE/IRS3C). 
IRS3E, which corresponds to #202, is classified as a low-mass CTTS, as also shown in figure 2 
in Preibisch et al. (2002). 

4-3. X-Ray Properties along the Evolutionary Phases 

We classified the 247 X-ray sources that have J-, H-, and A'-band counterparts (hereafter 
"NIR-bright sources") into CTTS and WTTS by the following criterion. Using the NIR data, 
we made the (J — H)/(H — K) color-color diagram shown in figure 6, where the colors were 
transformed into the CIT (California Institute of Technology) system (Bessel, Brett 1988) with 
a formula given by Carpenter (2001). The thin solid curves are the intrinsic colors of giant and 
dwarf stars (Tokunaga 2000), while the de- reddened locus (thick solid line) for CTTSs was taken 
from Meyer et al. (1997). We assumed the slope of the reddening line to be E(J — H)/E(H — 
K)= 1.63 (Martin, Whittet 1990), as given by the dashed lines. Since CTTSs are surrounded 
by cold accretion disks, they generally have larger A'-values than WTTSs. Therefore, we regard 
the sources located between the middle and the right reddening lines and above the CTTS locus 
as CTTSs, while sources located between the left and middle reddening lines are regarded as 
WTTSs. We obtained 57 and 164 sources with this classification. The remaining 26 sources 
located in other regions could not be classified. The expected errors in the classification are 
presented by Meyer et al. (1997). They concluded that approximately 1/3 of the CTTSs would 
occupy the region of the WTTSs. However, in this paper, we use the terms "CTTSs" and 
"WTTSs" for plainness. 

Excluding the unclassified sources, the number of bright sources and the averaged X-ray 
properties in each evolutionary phase are listed in table 3. We can see a higher kT for CTTSs 
(3.1 keV) than WTTSs (1.9 keV). However, there should be a selection bias in favor of a less- 
absorbed WTTS. Then, we also derived the averaged kT for the sources that had a Ah value 
less than 10 22 cm -2 . Nevertheless, the kT of CTTSs (2.4 keV) is still slightly higher than that 
of WTTSs (2.0 keV). The confidence level of the difference between the two distributions of kT 
is estimated to be 87.6% using the ASURV statistical software package (ver.1.2). 4 In order 
to examine this feature for other SFRs, we performed the same analysis using the published 
CXO catalogs of the Orion Nebula Cluster (hereafter ONC, Feigelson et al. 2002) and the 
Orion Molecular Cloud-2/3 (hereafter OMC-2/3, Tsujimoto et al. 2002), other sites of high- 
and intermediate-mass star formation. In identifying the 2MASS counterpart of X-ray sources 
in ONC, we defined the criterion that the angular separation of X-ray and NIR pairs is smaller 
than the solid line in figure 2, which is the envelope curve of the source identification criterion 
for Mon R2. For OMC-2/3, we used sources that were already identified with 2MASS sources 
by Tsujimoto et al. (2002). We then classified the sources with the same criterion as applied 
for Mon R2 and compared the mean X-ray temperature for each group, as shown in table 3. 

4 See (http://www.astro.psu.edu/statcodes/). 
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We again confirmed that the results similar to Mon R2; kT of CTTSs is slightly higher than 
that of WTTSs. The confidence levels of difference between the two distributions of kT are 
> 99.9% and 46.2%, respectively. Although the significance level is not high enough, except 
for ONC, we concluded that there is a marginal trend that kT of CTTSs is higher than that 
of WTTSs for all of the clouds (Mon R2, ONC, and OMC-2/3). This fact indicates that there 
is some changes between CTTSs and WTTSs in a physical process for the X-ray emission. In 
the case of YSO flares, Imanishi et al. (2003) reported that younger sources show a higher 
kT than evolved sources using a high-quality sample in the p Ophiuchi low-mass SFR. From 
a systematic study on the time profile of these flares, they concluded that the magnetic field 
gradually declines along the evolutionary phases, because the magnetic field strength has a 
positive correlation with kT (Shibata, Yokoyama 2002). For Mon R2, we can not discuss the 
flare and quiescent phases separately because of the limited photon counts. Future observations 
with a deep exposure time will reveal whether a similar effect can also be seen in Mon R2. 

4-4- X-ray Properties of Different Stellar Masses 

We divided the NIR-bright sources into two mass ranges, defined as M > 2.0 M & (high- 
and intermediate-mass), and M < 2.0 M & (low-mass) by the following criterion. In the J / (J — 
H) color-magnitude diagram (figure 7), the 1 Myr isochrone curve of 0.02 M Q to 3.0 M sources 
(Baraffe et al. 1998; Siess et al. 2000) is shown by the solid line. We assumed the slope of the 
reddening line to be Ajj E(J — H) =2.26 (Martin, Whittet 1990). We regard the sources located 
above the reddening line of stars with the mass of 2.0 M Q (dashed line) and rightward from 
the isochrone curve (solid curve) as high- and intermediate-mass sources, while sources below 
the reddening line and rightward from the isochrone curve to be low-mass sources. We then 
classified 10 and 230 sources into high- and intermediate-mass sources, and low-mass sources, 
respectively. The seven sources located leftward of the isochrone curve could not be classified. 
Excluding the seven unclassified sources, the number of bright sources and the averaged X-ray 
properties of each mass range are given in table 4. The Nn value for the low-mass sources are 
significantly smaller than that of the high- and intermediate-mass sources. This tendency may 
be due to the bias effect; the low-mass sources with large absorption come more easily below 
the detection limit compared to those of the high- and intermediate-mass sources, because the 
former has a systematically lower Lx than the latter (Kohno et al. 2002). 

Since we found a high kT (2.7 keV) of the high- and intermediate- mass sources similar 
to that of the low-mass sources (2.0 keV), we also performed the same analysis for ONC and 
OMC-2/3, and found similar results to Mon R2 (table 4). Moreover, a significant fraction of the 
high- and intermediate-mass sources show a rapid time variability. These results indicate that 
the X-ray emitting process of high- and intermediate-mass YSOs is the same as that of the low- 
mass sources, i.e., magnetic activity (Kohno et al. 2002; Preibisch et al. 2002). Because high- 
and intermediate-mass YSOs lack surface convection layers (Palla, Stahler 1993), the convective 
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dynamo process may not work. What kind of the magnetic amplification mechanism is working 
in high- and intermediate-mass YSOs is still an open question. 

4-5. The Spatial Distribution of X- Ray Sources 

The spatial distribution of X-ray sources in Mon R2 seems to be elongated in the same 
direction as the 13 CO intensity structure, from north to south (figure 1). We examined whether 
or not the X-ray and NIR source distributions have the same structures as that of 13 CO by 
evaluating the correlation coefficient R (Bevington, Robinson 1992). We defined two source 
groups: those detected by X-rays and NIR. We set the x- and y-axis along the edge of the 
ACIS-I FOV, and calculated R for each source group separately. We then obtained R = —0.19 
and —0.04 for the X-ray detected group and the NIR-detected group, respectively. In order 
to estimate the significance level, we randomly distributed the same number of test sources as 
those of X-ray and NIR in the x-y plane by 10000 times, and evaluated R. We then estimated 
the significance levels to be 3.7 a and 1.3 a, for the X-ray and NIR sources, respectively. We also 
derived R for sources that are detected by NIR, but not by X-ray, and obtained R = —0.03, the 
significance level of which is 1.0 a. Considering the significance levels, the distribution of the 
X-ray sources has a more elongated structure than that of NIR. This supports the observational 
review of Feigelson and Montmerle (1999) asserting that X-rays are a homogeneous tracer of 
all stages of YSOs, while the X-ray emission decreases as sources are approaching to the main- 
sequence stage. We confirmed this fact quantitatively for the first time by comparing the spatial 
distribution of X-ray sources and NIR sources. 

5. Summary 

We summarize the analysis of the 100 ks observation on Mon R2 with CXO ACIS-I as 
follows: 

1. We detected 368 X-ray sources from the cloud and performed spectral and timing analyses 
for most of the sources. This provided one of the largest catalogs of X-ray sources in an 
SFR. 

2. About 80% of the detected X-ray sources have NIR counterparts in the catalogs of 2MASS, 
SQIID, and IRCAM3. From their A% and i^-band magnitudes, we conclude that most of 
these NIR-identified X-ray sources are likely cloud members. 

3. We classified the cloud members into CTTSs and WTTSs using a (J — H)/(H — K) color- 
color diagram and compared the fractions of the X-ray time variability, averaged absorp- 
tion, and temperature. We found a marginal trend that the temperatures of CTTSs are 
slightly higher than those of WTTSs. We performed the same analysis for the ONC and 
OMC-2/3 members, and obtained similar results as Mon R2. Thus, the physical structure 
of X-ray emission may change from CTTS to WTTS. 

4. We divided the cloud members into high-, intermediate-, and low-mass sources based on 
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a J/(J — H) color-magnitude diagram, and compared the averaged X-ray properties. A 
significant fraction of the high- and intermediate-mass sources are time variable and exhibit 
a high-temperature plasma as those of low-mass sources. This supports the earlier results 
by Kohno et al. (2002) and Preibisch et al. (2002) that high- and intermediate-mass YSOs 
emit X-rays via magnetic activity. 
5. We performed a imaging analysis of the entire cloud core and the NIR cluster of Mon R2. 
We then found a significant difference in the spatial distribution between X-ray and NIR 
sources. We confirmed that X-rays can be used as a homogeneous tracer of all the stages 
of YSOs, as asserted by Feigelson and Montmerle (1999), and that the X-ray emission 
decreases as sources are approaching the main-sequence stage. 
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Fig. 1. ACIS-I image of the entire cloud core and the NIR cluster of Mon R2 (gray scale) overlaid on the 
13 CO ( J—l— 0) intensity contour map with a beam size of 2'6 in diameter (Loren 1977). The square shows 
the ACIST FOV of 17.'4 x YJ'A. The background photons are extracted from the circular region. Contours 
are the antenna temperature at the intervals of 1 K. The plus indicates the position of the maximum 
temperature. 



12 




2 4 6 8 10 12 



off-axis angle (arcmin) 

Fig. 2. Angular separation between X-ray sources and the NIR counterparts plotted as a function of the 
off-axis angle. The identification criterion for the ONC sources is shown by the solid line. 
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Fig. 3. Histogram of Nn of the bright sources. The mean value and that weighted with errors of the 
sources arc shown. 
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Fig. 4. Same as figure 3, but for kT. 
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Fig. 6. (J — H) I (H — K) color-color diagram in the CIT color system. X-ray sources that have J-, H-, 
and -RT-band counterparts are plotted with circles. The thick solid line indicates the de-reddened CTTS 
locus (Meyer et al. 1997). The intrinsic colors of giants and dwarfs are given by the thin solid curves 
(Tokunaga 2000). Reddening lines are given in dash (Martin, Whittct 1990). 



14 



-2 



CO 

E 



2 - 



~> 4 - 







o 
w 

.Q 

CD 



6 - 



8 - 



10 



3.0 Msun 




typical error 
Av=1 mag 



high- and 

intermediate-mass 



0.02 Msun 



low-mass 



J I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I L 



0.5 1 1.5 2 2.5 3 3.5 4 
J-H (mag) 



Fig. 7. Same as figure 6, but J /{J — H) color-magnitude diagram. The term "absolute J" indicates the 
J-band magnitude of each source is not corrected against the Ay value. The solid curve is the 1 Myr 
isochrone with 0.02 M Q < M < 3.0 M Q (Baraffe et al. 1998; Siess et al. 2000). The reddening line (dashed) 
is for stars with the mass of 2.0 M Q (Martin, Whittet 1990). 
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Table 1. Chandra X-ray sources and stellar counterparts. 
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1 Background-subtracted X-ray photon counts in the 0.5-10.0 keV band. 

* X-ray counterparts; (KKH): Kohno ct al. (2002) and (GMCF): Gregorio-Hetem et al. (1998). 
§ NIR counterparts; (CMDSH): Carpenter et al. (1997) and (2MASSI): 2MASS Point Source 
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Table 2. Results of timing and spectral analyses for bright sources. 
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Table 2. (Continued) 
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1.8 


(0.9-5.8) 


54.6 


^54.5-54.6) 


31.5 


161 


0.77 


0.2 


(0.0-1.0) 


1.3 


(0.5-1.9) 


52.6 


(52.5-52.7) 


29.6 


162 


0.60 


2.6 


(1.4-3.3) 


2.9 


(1.8-10.4) 


53.3 


(53.2-53.4) 


30.4 


165 


0.99 


8.9 


(5.6-15.6) 


2.0 


(0.9-3.6) 


54.2 


'54.1-54.2) 


31.2 


168 


1.00 


2.4 


(1.4-3.8) 


5.3 


(2.2-26.0) 


53.3 


(53.2-53.4) 


30.5 


169 


0.98 


3.9 


(2.5-6.0) 


2.5 


(1.2-...) 


53.7 


(53.6-53.8) 


30.7 


172 


0.48 


0.8 


(0.2-1.4) 


1.4 


(0.9-3.3) 


52.9 


(52.8-53.0) 


29.9 


179 


0.99 


6.5 


(3.6-11.3) 


2.6 


(1.2-8.3) 


53.7 


(53.6-53.8) 


30.8 


183 


0.88 


1.5 


(0.9-2.6) 


2.6 


(1.5-6.9) 


53.1 


(53.0-53.2) 


30.1 


184 


0.94 


0.7 


(0.2-1.7) 


2.0 


(1.0-5.7) 


52.8 


(52.7-52.9) 


29.8 


190 


0.99 


0.6 


(0.3-1.4) 


13.0 


(...-...) 


52.7 


(52.6-52.8) 


30.0 


192 


0.99 


0.4 


(0.2-0.8) 


3.0 


(1.7-...) 


53.2 


(53.2-53.3) 


30.3 


193 


0.60 


4.8 


(3.0-7.0) 


2.1 


(1.1-4.3) 


53.9 


(53.8-53.9) 


30.9 


197 


1.00 


10.5 


(5.6-20.1) 


2.0 


(...-8.3) 


54.0 


(53.9-54.1) 


31.0 


198 


1.00 


1.7 


(1.2-2.5) 


5.1 


(2.8-...) 


53.6 


(53.5-53.6) 


30.7 


202 


0.99 


0.8 


(0.2-2.4) 


3.9 


(...-...) 


52.8 


(52.7-52.9) 


30.0 


209 


1.00 


2.7 


(2.5-2.9) 


4.3 


(3.6-5.2) 


54.7 


(54.7-54.7) 


31.8 


210 


0.99 


1.0 


(...-2.2) 


3.4 


(...-62.0) 


53.3 


(53.2-53.4) 


30.4 


212 


0.97 


5.0 


(3.2-...) 


2.8 


(...-...) 


53.9 


(53.9-54.0) 


31.0 


213 


0.33 


1.0 


(0.7-1.3) 


2.0 


(...-...) 


53.5 


(53.5-53.6) 


30.5 


217 


0.99 


0.3 


(0.2-0.4) 


1.3 


(1.1-1.5) 


53.4 


(53.3-53.4) 


30.3 


220 


0.98 


0.3 


(0.1-0.6) 


1.9 


(1.3-2.8) 


53.0 


(52.9-53.0) 


29.9 


224 


1.00 


2.4 


(1.4-3.6) 


2.9 


/in \ 

(1.6-...) 


53.6 


(53.5-53.7) 


30.7 


225 


0.93 


1.8 


(1.2-2.8) 


3.6 


(1.9-10.7) 


53.2 


(53.1-53.3) 


30.3 


226 


0.95 


1.2 


(0.4-1.8) 


1.1 


(0.7-1.6) 


53.3 


(53.2-53.4) 


30.2 


228 


0.99 


1.3 


(0.8-1.8) 


1.1 


(0.8-1.5) 


53.5 


(53.4-53.5) 


30.4 


229 


0.97 


0.7 


(0.4-1.3) 


4.1 


(...-15.5) 


53.1 


(53.0-53.2) 


30.3 


231 


0.85 


0.8 


(...-...) 


0.9 


(...-...) 


53.4 


(53.3-53.5) 


30.4 


237 


0.91 


1.8 


(0.7-2.6) 


1.2 


(0.7-1.7) 


53.4 


(53.3-53.5) 


30.3 


255 


0.85 


1.0 


(0.5-1.9) 


2.2 


(1.1-4.6) 


53.0 


(52.9-53.1) 


30.0 


256 


0.79 


2.2 


(1.5-3.5) 


10.0 




53.1 


(53.0-53.2) 


30.2 


259 


0.92 


0.8 


(0.2-...) 


0.9 


(...-...) 


53.5 


(53.4-53.6) 


30.5 


261 


0.72 


2.2 


(1.6-3.6) 


1.7 


(0.9-3.2) 


53.9 


(53.8-54.0) 


30.9 


264 


0.98 


1.2 


(0.5-2.6) 


4.0 


(...-60.8) 


53.0 


(52.9-53.1) 


30.1 


265 


0.43 


0.7 


(0.5-1.0) 


0.7 


(...-...) 


53.5 


(53.4-53.5) 


30.4 


268 


0.71 


1.2 


(0.5-1.6) 


0.4 


(...-2.4) 


53.7 


(53.6-53.8) 


30.6 


270 


0.46 


0.9 


(0.1-1.4) 


1.0 


(0.5-11.0) 


53.0 


(52.9-53.1) 


30.0 


277 


0.99 


2.9 


(...-...) 


1.7 


(...-...) 


53.7 


(53.7-53.8) 


30.7 


Zoi 


U.yo 


0.1 


(0.1-0.4) 


A Q 
4.0 


(A A \ 

(1.4-...) 




oz.o— oz.o ) 


Zy.o 


286 


0.97 


1.2 


(0.8-2.1) 


2.9 


(1.4-7.3) 


53.2 


(53.1-53.3) 


30.2 


289 


0.94 


1.2 


(...-1.4) 


1.1 


(...-...) 


54.2 


(54.2-54.2) 


31.1 


293 


0.95 


0.4 


(0.2-0.9) 


4.2 


(2.0-12.3) 


52.8 


(52.7-52.9) 


30.0 


295 


0.84 


2.2 


(1.5-3.3) 


10.0 




53.1 


(52.9-53.2) 


30.2 


296 


1.00 


0.3 


(0.1-0.4) 


2.3 


(1.6-3.3) 


53.2 


(53.2-53.3) 


30.2 


301 


0.97 


0.5 


(0.1-...) 


0.7 


(0.5-1.1) 


53.5 


(53.5-53.6) 


30.5 


304 


0.99 


0.6 


(0.4-...) 


2.9 


(2.3-3.4) 


53.9 


(53.9-53.9) 


31.0 
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Table 2. (Continued) 



ID Variable 
probability* 



(10 22 cm 



2 ) 



keV) 



3 ) 



cm 



logOLx 11 ) 1 
(erg s" 1 ) 



305 
307 
309 
310 
311 
313 
317 
320 
321 
329 
331 
333 
344 
347 
349 
350 
351 
353 
354 
355 
357 
358 
359 
360 
363 
366 
367 



0.94 
0.82 
0.99 
0.99 
0.99 
0.94 
0.99 
0.99 
0.92 
0.78 
0.82 
0.99 
0.96 
1.00 
0.99 
1.00 
0.99 
0.60 
0.99 
0.56 
0.88 
0.99 
0.96 
1.00 
0.99 
0.63 
0.91 



1.4 
0.6 
0.3 
0.5 
0.6 
0.2 
0.6 
0.2 
1.4 
0.4 
3.6 
0.3 
0.2 
0.2 
1.7 
1.3 
0.6 
1.0 
1.0 
0.8 
4.2 
0.8 
1.0 
0.5 
2.2 
0.8 
0.6 



1. 
0. 
0.1 



0- 1.9 

1- 1.1 
0.7 

•) 



-0.4 

-■■•) 
-0.4 
-2.1 
-0.7 





2.6-5.5 



-0.8 
-0.9 
-0.3 
-2.8 
-■■•) 



1 
1 
1 
1 
2 

2-1.1 
7-1.5 
2-1.6 

5- 1.1 

6- 6.6 
-...) 

2.0 
0.9 
2.5 
...) 
1.4) 



20.2 
0.9 
3.4 
1.6 
1.5 
1.6 
1.9 
2.1 
4.9 
1.3 

23.1 
1.5 
1.2 

0.9 
4.3 
0.9 

10.0 
1.0 
0.5 

10.0 
1.0 

28.3 
3.7 
2.8 
1.7 
2.1 



...-...) 

0.6-1.7) 

1.7-8.3) 

1.3-1.8) 

1.3-1.8) 

1.0- 2.9) 

1.1- 2.4) 
1.6-...) 
...-11. 5) 
1.0-1.7) 
...-...) 
0.8-2.5) 
0.5-1.6) 

0.5-1.9) 
3.6-5.2) 
0.7-1.4) 

0.7-1.7) 
...-...) 

0.8-2.3) 
...-...) 
...-22.7) 
...-3.5) 
1.3-2.0) 
1.0-3.1) 



53.8 
53.0 
52.8 
54.2 
54.0 
53.2 
53.4 
53.5 
53.4 
53.1 
53.7 
53.2 
52.8 

53.4 
54.5 
53.1 
53.2 
53.2 
53.3 
53.2 
53.7 
52.9 
52.8 
54.4 
54.1 
53.2 



53.8- 53.8 

52.9- 53.1 



52.7- 
54.1- 



52.8 
54.2 



54.0- 54.1 

53.1- 53.3 

53.4- 53.5 

53.5- 53.6 

53.4- 53.5 

53.0- 53.1 

53.6- 53.7 

53.1- 53.3 

52.7- 52.9 

53.3-53.5 

54.5- 54.6 

53.0- 53.2 

53.1- 53.3 



53.1- 
53.1- 



53.3 
53.4 



53.0- 53.4 

53.7- 53.8 

52.8- 53.0 
52.7-52.9 
54.3-54.4 

54.1- 54.1 
53.1-53.3 



31.1 
30.0 
29.9 
31.1 
31.0 
30.2 
30.4 
30.5 
30.6 
30.0 
31.0 
30.2 
29.8 
29.8 
30.4 
31.7 
30.0 
30.3 
30.1 
30.2 
30.3 
30.7 
30.2 
30.0 , 
31.4(31 
31.0 (31 
30.2 (30 



(31. 
(29. 
(29. 
(31 
(31 
(30. 
(30. 
(30. 
(30. 
(29. 
(30. 
(30. 
(29. 
(29. 
(30. 
(31 
(29. 
(30. 
(30. 
(30. 
(30. 
(30. 
(30. 
(29 



0- 31.1) 
9-30.1) 
8-29.9) 

1- 31.1) 
0-31.0) 



30.3) 
30.4) 
30.6) 
30.7) 
9-30.1) 
9-31.0) 

1- 30.2) 
7-29. 
7-29.8) 

2- 30.5) 
7-31.7) 
9-30.1) 
2-30.4) 
0-30.2) 



1-30.3) 
1-30.5) 
6-30.7) 
0-30.2) 
9-30.1) 
4-31.4) 

0- 31.1) 

1- 30.3) 



~* Variable probabilities obtained by a Kolmogorov-Smirnov test to the photons in the 
0.5-10.0 keV band. 

T 90% confidence limits are in parentheses. No value in a parenthesis means that the error 
is not constrained. 

* If sources have no converged value of kT, we assume that they have fcT=10.0 keV. For 
#347, the thin-thermal plasma model is rejected, while power-law model is accepted. 

§ Emission measures. 

II Luminosity in the 0.5—10.0 keV band. 



19 



Table 3. X-ray properties along the evolutionary phases. 





r lr r r rQc! 
iss 


W 1 1 OS 


JVLOII Kz 






JN umber 


z3 (12) 


69 (45) 


Variability^ (ratio %) 


O f O A 0\ 


23 (33.3) 


at hn22 ni „-2U 
-^*H cm J + 


i q /i 
l.o t^l-UJ 


u.y (^u.oj 


Kl i yKQ V J + 


o.l \^-o) 


i.y \lA) 




2.4 (1.1) 


z.U (1.5) 








Number 1 * 1 


DO (00J 


zoo ^zOU J 


Kl i t^xe V J + 


z.y (i.oj 


O Q /I /j \ 
Z.O ^1.4J 


kT 2 (keV)+ 


3.1 (1.6) 


2.4 (1.4) 


OMC-2/3 






Number* 


14 (7) 


92 (70) 


kTi (keV)* 


2.1 (1.1) 


1.4 (0.7) 


kT 2 (keV)i 


1.5 (0.7) 


1.3 (0.7) 



* The number of sources that have iVjj value of less than 10 



cm -2 are given in parentheses. In the case of Mon R2, the 
number of the bright sources in each evolutionary phase are 
shown. 

t The number of the time-variable sources. The fraction of the 
number of time- variable sources are given in parentheses, 
t The average of the best fit values (not weighted by the errors) . 
The standard deviations are given in parentheses. 
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Table 4. X-ray properties of different stellar masses. 

High- and Low mass 
intermediate- 



Mon R2 
Number* 

Variability^ (ratio %) 
N H (10 22 c 



-2U 



kT (keV)l 



3 (37.5) 

2.6 (2.6) 

2.7 (1.9) 



92 

33 (35.9) 
1.0 (0.8) 
2.0 (1.4) 



ONC 
Number 
kT (keV)* 



14 

3.1 (1.4) 



372 

2.4 (1.4) 



OMC-2/3 

Number 9 94 

kT (keV)* 2.1 (1.2) 1.4 (1.2) 

* The number of the bright sources in each mass range. 

t The number of the time-variable sources. The fraction of the 

number of time- variable sources are given in parentheses. 

t The average of the best-fit values (not weighted by the errors). 

The standard deviations are given in parentheses. 
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